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We report on a molecular dynamics investigation of tlie wetting properties of graphitic surfaces by various 
solutions at concentrations 1 — 8 wt% of commercially available non-ionic surfactants with long hydrophilic 
chains, linear or T-shaped. These are surfactants of length up to 160 [A]. It turns out that molecular dynamics 
simulations of such systems ask for a number of solvent particles that can be reached without seriously com- 
promising computational efficiency only by employing a coarse-grained model. The MARTINI force field with 
polarizable water offers a framework particularly suited for the parameterization of our systems. In general, 
its advantages over other coarse-grained models are the possibility to explore faster long time scales and the 
wider range of applicability. Although the accuracy is sometimes put under question, the results for the wetting 
properties by pure water are in good agreement with those for the corresponding atomistic systems and theo- 
retical predictions. On the other hand, the bulk properties of various aqueous surfactant solutions indicate that 
the micellar formation process is too strong. For this reason, a typical experimental configuration is better ap- 
proached by preparing the droplets with the surfactants arranged in the initial state in the vicinity of contact line. 
Cross-comparisons are possible and illuminating, but equilibrium contanct angles as obtained from simulations 
overestimate the experimental results. Nevertheless, our findings can provide guidelines for the preliminary 
assessment and screening of surfactants. Most importantly, it is found that the wetting properties mainly depend 
on the length and apolarity of the hydrophobic tail, for linear surfactants, and the length of the hydrophilic head- 
group for T-shaped surfactants. Moreover, the T-shaped topology appears to favor the adsorption of surfactants 
onto the graphitic surface and faster spreading. 

PACS numbers: 47.11.Mn,47.55.dk,47.55.dr,47.55.np 



I. INTRODUCTION 



Molecular dynamics (MD) simulations have become an 
important complement to experimental and theoretical work 
across a variety of disciplines. Often, one of the encoun- 
tered obstacles to further advances are the prohibitive compu- 
tational resources necessary for a description at the molecular 
level of systems containing a large number of particles |1J. 
This situation is in general quite common when water is in- 
volved as solvent. One way to overcome this limitation is to 
reduce the degrees of freedom to the essential for the issues 
under investigation. This procedure is referred to as coarse- 
graining because the system is suitably organized into units 
identifying groups of atoms. To achieve this, two schemes 
are basically available |2--^. The differences mainly reside 
in the decomposition of the system into building blocks and 
the selected properties according to which inter-particle inter- 
actions are calibrated. Our work is based on the MARTINI 
force field because allowing to build readily and systemat- 
ically systems of coarse-grained (CG) units in a prescribed 
way |2 |. Within this framework four heavy atoms are gener- 
ally mapped into one interaction site. The resulting beads are 
classified into four main categories that separate further into 
subtypes. Lennard-Jones (LJ) interactions are organized into 
ten levels of strength and evaluated with a function of type 
12-6. In its original formulation, the MARTINI model does 
not contemplate solid phases |2|. Following empiiical argu- 
ments, we treat graphitic surfaces and their wetting properties 
by surfactant solutions. In that respect, the surface tension 
of water is of course a key property. The polarizable MAR- 
TINI model for water |3 1 still underestimates the experimental 



value. The expected balance of forces at the interface with the 
graphitic surface is approximately restored by tuning the in- 
teraction strength of the relative LJ potential. In this way, we 
first reproduce in the macroscopic limit the contact angle for 
pure water from measurements at present carried out under the 
most ideal conditions [6j. Interestingly, our results correlate 
well with those using atomistic models for water consistent 
with the experimental surface tension 111 [SI. This allows us 
to conclude that the MARTINI force field can provide an ef- 
fective framework for investigating the wetting properties of 
surfaces. Surfactants are molecules exhibiting opposite be- 
havior toward water at their extremes, namely hydrophilic in 
the headgroup and hydrophobic in the tail group. Importantly, 
these competing attributes are responsible for the decrease of 
the surface tension of the resulting solution 19]. In practice, 
this means that an aqueous solution would wet better a surface 
with the addition of surfactants. At higher concentrations, an- 
other important property related to the amphiphilic nature of 
surfactants is their ability to self-assemble into micelles |9|. 
Applications in the industry are innumerable, including print- 
ing, detergency, oil recovery, dispersion, emulsification, dis- 
infection, etc. H. In this work, we concentrate on solutions 
with non-ionic surfactants. As conveyed by the title, empha- 
sis is put on the wetting propeities by measuring the contact 
angle of droplets. The identity of the surfactants under investi- 
gation is confidential for industrial reasons. Specifically, their 
use as wetting agents is necessary for the optimal dispersion 
of reinforcement fillers in polymeric matrices (low chemical 
compatibility to each other) This phase is rather promi- 
ment in the manufacturing route of composite materials in or- 
der to meet challenging, target properties llTTIfTSl . Our work 
aims at providing useful guidance for planning applications. 
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So, the basic mechanisms governing wetting phenomena by 
aqueous surfactant solutions are at the center of our attention. 
The strategy consists in simulating six representative surfac- 
tants differing in simple attributes such as the length of the 
hydrophilic and hydrophobic groups or the topology (linear 
and T-shaped). Our research develops as follows. In the next 
Section, the MARTINI representation of all the constituents 



of the systems is introduced. Section III is concerned with the 



general MD simulation settings. In Sec. IV are presented the 
results for systems of pure water and with surfactants: first the 
bulk properties of solutions and then the wetting behavior on 
graphitic surfaces. The last Section is devoted to conclusive 
remarks. 



II. COARSE-GRAINING OF THE SYSTEMS 

Force field. The MARTINI force field is designed in order 
to account for a variety of chemical compounds preventing the 
need to recalibrate the various interaction parameters for ev- 
ery application f2}^. The idea is to classify atomic functional 
groups into four main categories: polar, nonpolar, apolar and 
charged. Every category subdivides further into subtypes so 
as to represent more finely the underlying structural proper- 
ties. Ten levels of interaction are possible between the diverse 
types of CG particles. More precisely, inter-particle LJ forces 
are derived from the potential energy 



t/Lj('-//) =4e,- 



The indices / and j label the particles. r,y is the distance sepa- 
rating pairs of particles. The parameter £,y adjusts the strength 
of the forces according to the available levels of interaction. 
The coefficient Oij is related to the equilibrium distance be- 
tween particles via the formula 2^6 (j,.^.. If • < l^l^Oij the 
force is repulsive, otherwise it is attractive. The parameter Oij 
is also usually referred to as molecular diameter Electrostatic 
interactions are computed from the coulombic potential 



qiqi 



47reoerr; 



qi is the partial charge of the /-th particle; Eq is the permit- 
tivity in vacuum and £r is the relative dielectric constant. For 
the standard MARTINI model = 15, while for the polaiiz- 
able model used here = 2.5. Non-bonded interactions are 
corrected separately to zero using the shift function lfT3l 



nj < '•shift 



3 inj - '"shift) + 4 ('",■,/■ - rsliift) + C rshift < nj < rc , 

where rghift is the inner cutoff distance and is the outer cut- 
off distance. The constants A, B and C are determined from 
boundary conditions. Intra-particle interactions take into ac- 
count bond and angle interactions. Bonded particles interact 
via the harmonic potential 



where ro is the equilibrium distance and K\, the energy con- 
stant. The standard MARTINI values for these parameters are 
fQ = 4.7 [A] and K\, = 1.5 [kcal/mol/A^]. Angle interactions 
are described by the potential 

=^a[cos(0„-i)-cos(0o)]^ 

Ojjk is the angle formed by triplets of bonded particles; Oq is 
the equilibrium angle; is the coupling constant, for which 
the standard MARTINI value is K^ = 3.0 [kcal/mol]. 

Water. A CG bead of polai'izable MARTINI water |3 1 con- 
sists of three particles Wq, W± of equal mass 24 [g/mol] and 
of partial charge qo — Q [e] and q± = ±0.46 [e], respec- 
tively. The charged particles are bonded to the neutral one 
and the bond length is kept fixed at 1.4 [A]. The angle de- 
gree of freedom is modeled as a harmonic oscillator of elas- 
tic constant = 0.5019 [kcal/mol] and equilibrium angle 
of 0°. Lennard-Jones forces are considered only among the 
central particles Wo with interaction parameters Cww = 0.956 
[kcal/mol] and Oww ~ 4.7 [A]. With the charged particles W±, 
the polarization is treated explicitly and the relative dielectric 
constant is thus set to Ci- = 2.5 (cf. Refs. ElO). Furthermore, 
Coulomb and van der Waals forces among particles belonging 
to the same water bead are omitted. 

Graphene. The graphene is parameterized by placing at the 
center of every ring of carbon atoms one interaction site. As a 
consequence, every bead has mass 24 [g/mol] because repre- 
senting on average two atoms. The interaction parameter ecw 
between the beads of water (Wo) and graphene is fixed from 
the wetting properties of pure water. We choose the value 
of Ecw that reproduces the result for the contact angle in the 
macroscopic limit reported in Ref. |6|. The parameter Ocw is 
set equal to Ocw — 6.24 [A]. This value yields approximately 
the equilibrium distance of 7 [A]. The solid phase is obtained 
by two parallel planes of CG graphene. The two planes are 
separated by 3.4 [A] with the lower one translated by the vec- 
tor {1/2, \/3//2) with respect to the upper one. Here / is the 
length of C-C bonds in the all-atom case. This arrangement is 
intended to reproduce the structure of a graphite crystal. 

Surfactants. The atomic structure of the surfactants is re- 
duced according to the guidelines detailed in Ref. Q . In order 
to compare our work to previous similar studies [T^l, we also 
consider solutions containing the non-ionic surfactant C8E4 
(tetraethylene glycol octyl ether). We adopt the mappings Ci- 
C1-P5-P5 and Ci-Ci-Na-Na-P4 proposed in Ref [14]. The 
first representation will be referenced as C8E4-P5 surfactant 
while the latter as C8E4-NaNaP4 1 14|. Both molecules are lin- 
ear with every equilibrium bond length equal to ro — 4.7 [A]; 
all equilibrium angles are of course 0o — 180°. The energy 
constants are also the standard MARTINI values [2], that is, 
K^ = l.5 [kcal/mol/A2] and/Ta = 3.0 [kcal/mol]. For the mass 
of the beads, we use in this case mp^ = 72 [g/mol], mNa ~ 56 
[g/mol] and mc, = 56 [g/mol]. The LJ parameters are read- 
ily obtained from Ref |2 1 with all equilibrium distances set to 
4.7 [A]. The wetting properties of graphitic surfaces are stud- 
ied in relation to five non-ionic surfactants, long-chain and 
T-shaped. The first surfactant will be referred to as LI. Its 
molecule is linear and it is applied the MARTINI representa- 
tion (Ci)3(P4)io. The mass is again mc; = 56 [g/mol] for the 
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tail beads, while we use mp^ = 44 [g/mol] for the hydrophilic 
head. The second surfactant has representation (Ci)3(P4)20 
and will be designated by L2. All topology and interaction 
parameters are the same as for the surfactant LI. The third 
linear surfactant L3 is mapped into (Ci)4(P4)3o. All the pa- 
rameters are the same as for the surfactant LI. Regarding the 
T-shaped surfactants, we distinguish the Tl and T2 topolo- 
gies. These surfactants differ by the length of the hydrophilic 
headgroup. The surfactant Tl has a linear chain made of ten 
beads P4, while the surfactant T2 counts twenty beads P4 ar- 
ranged linearly. In both cases, the hydrophobic tail consists of 
three aligned beads Ci and in the middle is attached the hy- 
drophilic chain. For the angle C1-C1-P4 (vertex in the central 
bead), the equilibrium value is set to 90°, and the standard po- 
tential and interaction parameter of the MARTINI force field 
1^ are used. All the other parameters are the same as for the 
linear surfactants. The hydrophobic tail of the surfactant T3 is 
made of five Ci beads. At the central bead is attached the lin- 
ear, hydrophilic head of length twenty P4 units. All the other 
settings for T-shaped surfactants apply also in this case. 



m. SIMULATIONS 

All simulations are performed with the molecular dynam- 
ics code LAMMPS ifTsHTTll . Newton's equations of motion 
are integrated with the Verlet algorithm using a Nose-Hoover 
scheme in the specified ensemble; the timestep size is of 20 
[fs]. Non-bonded interactions are cut off at rc = 12 [A] and 
computed with the potential Ij/gromacs/coul/gromacs lfT3l . 
The LJ potential is shifted from rjhift = 9 [A] to r^, while the 
electrostatic contribution is shifted from rshift = [A] to rc. 
Long-range interactions are not taken into account. The bond 
length between Wo and W± particles is maintained rigid with 
the SHAKE algorithm |18|. Three -body interactions of sur- 
factants are computed with the cosine/squared potential. The 
neighbor list is rebuilt at most every 5 timesteps. It is not nec- 
essary to divide the bond and angle coefficients given here by 
2. Then, all simulation times must be meant as actual |3 |. Un- 
less specified otherwise, these parameters remain unchanged 
in the various simulations. The motivations for these settings 
can be found in Refs. 12) (3] • 



IV. RESULTS AND DISCUSSION 

Bulk properties of CG water. As starting configuration, 
I'OOO beads of water are arranged with the Wo particles on 
the vertices of a simple cubic lattice. The side of the unit cell 
is 4.5 [A]. The system is let evolve for 100 [ns] at NPT condi- 
tions. The time of evolution is long in order to verify that we 
do not incur in freezing. The target temperature and pressure 
are r = 298 [K] and P ~ \ [atm]. The system is studied by 
recording I'OOO evenly-spaced frames over the course of the 
last 50 [ns]. Table |l] lists some average bulk properties of the 
cubic box of CG water The Wo-Wo radial distribution func- 



L[A] 


V [A3] 


P [g/A^] 


d [molecules/A-'] 


48.858 


11 6' 629 


L0255 


0.0086 



Table I: Average characteristics of the cubic simulation box during 
the NPT dynamics, side length L and volume V, along with the bulk 
properties of CG water p and d, mass and particle densities respec- 
tively. 
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Figure 1: Wq-Wo RDF as obtained from the NPT dynamics for the 
polarizable MARTINI water model. 



tion, g(r), is computed by means of the formula 

^g{r)A%r^Ar = S{r) . 

N is the number of water beads and V the volume of the sim- 
ulation domain. S{r) counts the average number of Wo parti- 
cles falling in a shell of width Ar = 0.05 [A] centered around a 
given Wo particle at distance r. The radial distribution function 
(RDF) is shown in Fig.[T] In general, our results compare well 
with those of the original work for the present model of CG 
water O. The standard deviation of the volume of the simu- 
lation domain is 0.78% of the average value. For this reason, 
we safely replicate the last configuration of this NPT simu- 
lation and extract the asseblies for any subsequent dynamics 
containing pure CG water. 

Spherical droplets. As a matter of calibration of the in- 
teraction parameter Ecw, droplets of pure CG water spread- 
ing on the CG graphene are first considered. The interaction 
strength ecw is varied in the range 0.04 — 0.08 [kcal/mol]. 
Every initial configuration is composed by a hemispherical 
droplet of 17'927 beads centered above the upper-lying plane 
of graphene at a distance of 7 [A]. The radius of the hemi- 
sphere is 100 [A]. The planes of graphene are squared with 
side length 320 [A]. The boundary conditions are periodic. 
Every system is evolved in the canonical ensemble for 60 [ns] 
at r = 298 [K]. The systems are studied from I'OOO frames 
recorded during the last 40 [ns]. These dynamics are shorter 
because for long times of evolution the droplets systematically 
exhibit excessive stratification. The onset of this phenomenon 
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Figure 2: Contact angle as a function of the strength parameter. 
Matching with the experimental result 16J occurs approximately for 
ecw = 0.0575 [kca^mol]. 
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Figure 3: (Color online) Profile of the spherical droplet for ecw = 
0.0575 [kcal/mol]. The data are analyzed following the method pre- 
sented in Refs. f8l[l9l|20|. The green line is a fit to the data under 
the assumption of a spherical shape of the droplet. The tangent at the 
contact line, where the solid, liquid and vapor phases meet, is in red 
and the contact angle is measured to be 128.2°. 



depends on several parameters, first of all the equilibrium dis- 
tance of the liquid interface from graphene. Presumably, this 
shortcoming might be avoided with an involved work of re- 
finement of the force field. We do not insist on this aspect 
because we expect that the presence of surfactant molecules 
would reduce the formation of layers of solvent particles. Fur- 
thermore, higher temperatures could retard this phenomenon 
because of faster kinetics. Then, the behavior at room temper- 
ature could be extrapolated. For all simulations presented in 
this Article, stratification never occured and, as we shall see, 
our outcomes are still in reasonable agreement with the atom- 
istic results and the general theoretical predictions. Figure [2] 
shows the dependence of the contact angle on the parame- 
ter Ecw- These data clearly satisfy a linear relation as in the 
atomistic case |8|. For ecw — 0.0575 [kcal/mol], it is found 



r[A] 


no. water beads 


e[°] 


80 


25'933 


128.0° 


90 


32'742 


128.7° 


100 


40'430 


129.0° 


110 


48' 800 


129.9° 


120 


58'173 


130.4° 


130 


68'255 


130.4° 



Table II: Contact angle as a function of the initial radius of cylinders. 
The variations of 6 indicate that the macroscopic regime is reached 
for radii ~ 130 [A]. 



r[K] 


d [molecules/ A^] 


Tf [A] 




298 


0.0076 


100.3 


130.4° 


305 


0.0075 


101.8 


129.3° 


310 


0.0075 


101.1 


130.6° 


315 


0.0075 


101.3 


130.8° 


320 


0.0074 


101.4 


131.4° 


350 


0.0071 


102.0 


132.9° 


400 


0.0065 


104.1 


136.2° 



Table III: Variations from temperature rise for the cylinder of initial 
radius 130 [A]: particle density d, final radius and contact angle Q. 



a contact angle of 128.2°. This result reproduces the experi- 
mental value reported in Ref. ||6] . The profile of the droplet in 
this case is shown in Fig.|3] In principle, other coarse-grained 
force fields are available for such investigations W^S^- For ex- 
ample, in Ref. EH the wetting properties of pure water on a 
flat surface were investigated by taking a LJ function of type 
6-3 for the interaction between the liquid and solid phases. 
Using the same CG water, preliminary results indicate that 
the process of stratification on a molecular graphitic substrate 
occurs quite fast with the liquid-solid interactions described 
by a LJ potential of type 12-6. Ultimately, the assessment of 
other coarse-grained force fields for wetting studies warrants 
a detailed discussion, which is beyond the scope of the present 
study. 

Cylindrical droplets. We now want to investigate the depen- 
dence of the contact angle on size and temperature. To achieve 
this, six hemicylindrical droplets of radii from 80 to 130 [A] 
are considered. The x side of the simulation domain is always 
300 [A], corresponding to the height of every hemicylinder. 
The y side is always 120 [A] larger than the diameter of the 
droplets. In the initial configuration, the droplets are centered 
in the xy plane. All other simulation settings are the same em- 
ployed for spherical droplets. In Tab. [H] are listed the contact 
angles measured from the profiles of cylinders with increasing 
radii. The contact angles vary only of a few degrees and thus 
we conclude that size effects are negligible with radii ^130 
[A]. It is interesting to remark that the contact angle for cylin- 
ders is slightly larger than that for the spherical droplet, as 
observed for atomistic systems approaching the macroscopic 
regime |8|. For the largest cylinder, different temperatures 



are considered. The results of Tab. Ill suggests that variations 
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over a range of 100 [K] have no drastic effect. The resuhs 
reported here for the temperature dependence of contact an- 
gle are in line with those for atomistic systems |7|. In the 
following, cylindrical droplets are preferred to spherical ones 
because the simulations turn out to be accelerated |22|. More 
precisely, the CPU timings for the cylindrical droplet of ra- 
dius 80 [A] are comparable with those of spherical droplets, 
of radius 100 [A]. The former system contains 25'933 water 
beads while the latter 17'927. The reason for the speed-up 
is that parallelization based on spatial decomposition |17| is 
more efficient for cylinders because the projection of the par- 
ticles on the xy plane of the simulation domain gives a better 
coverage. 

Bulk properties of surfactant solutions. Here we consider 
boxes of CG water containing surfactants. In the initial con- 
figuration, the surfactant molecules are arranged regularly be- 
tween two slabs of water Every simulation lasted for 300 [ns] 
at NPT conditions. The target pressure is always 1 [atm] while 
the temperature is adjusted to different values. Data for anal- 
ysis are recorded every 5 [ns]. To start with, we want to com- 
pare the predictions of the MARTINI model with polarizable 
water |3| to those of the standard version |2|. To this end, we 
first simulate 175 surfactant molecules of C8E4-NaNaP4 with 
18'339 water beads, as done in Ref. |14|. This concentration 
leads to the formation of micelles lil4J . We look at the RDF 
between Ci beads belonging to different surfactant molecules 
in order to extract the equilibrium distance between tail beads 
in the same micelle. From Fig. |4] we see that there is a first 
peak around 5.1 [A], while the first minimum occurs at about 
'"mill = 7 [A]. As a consequence, we shall assume that the sur- 
factants with all of their C i beads separated by more than r^^m 
are free with respect to each other. The first pronounced peak, 
together with the fact that the RDF tends to zero, indicate that 
the tail beads are confined in relatively small regions of the 
whole simulation domain. Indeed, for all temperatures it is 
found that most of the surfactant molecules cluster into mi- 
celles, corroborating the results reported in Ref. |14|. Table 
HV] lists basic statistical indicators relative to the size distri- 
bution of micelles. As the temperature increases, the average 
size of micelles increases and there is clearly more disper- 
sion around the mean value (enhanced polydispersity). Fur- 
thermore, the size of micelles (aggregation number) tends to 
increase with temperature, as expected for non-ionic surfac- 
tants |23|. In that respect, it is interesting to remark that the 
first peak of the Ci-Ci (tail beads) RDF is more marked at 
higher temperatures (see Fig.|4]). For the other representation 
of the surfactant C8E4 1 14| (see Sec.[ll]i, we also consider 175 
molecules with 18'339 water beads at the temperature of 350 
[K]. The Ci-Ci RDF for the surfactant C8E4-P5 is similar to 
those of Fig. |4] and the characteristic distance r,mn associated 
with micelle formation is 7 [A] also in this case. For this map- 
ping, the headgroup is more polar and the process of micelle 
formation is effectively weakened. Figure |5] shows the num- 
ber of free surfactants in the course of time. Table IV presents 



micelles leaves the surfactant solution at CMC. After 300 [ns] 
of dynamics th ere remain five free surfactant molecules (see 
Fig.|5]and Tab.|iv]i: it is thus predicted a CMC of 4.4 x 10"^ 
[M]. Now, the CMCs at two different temperatures Ti and T2 
are related by CMC(r2) = [CMC(ri)]^i/^2 04.21]. From 
the experimental value CMC(ri = 298[K]) = 8 x 10"^ [M] 
lHHHg it is obtained a CMC of 16 x lO^^ at 350 [K]. It turns 
out that our result underestimates this value but is closer than 
the finding within the standard MARTINI model of 35 x 10"^ 
[M] [14]. For the surfactant L2, L3, Tl and T2, we apply the 
same simulation settings with the exception that the timestep 
size is reduced to 10 [fs] for stability reasons related to bond 
interactions of the surfactant molecules. The temperature is 
always kept fixed atT — 298 [K] and the initial state consists 
of a given number of surfactant molecules, for the concentra- 
tions of Tab. |v] between two slabs of water beads (50'000 in 
total). For the typical experimental concentration of 1 wt%, 
the RDFs of Fig. [6] indicate that the tail beads, for the sur- 
factant L2, are separated by larger distances because the first 
maximum is lower. From the size distributions of micelles 
plotted in Fig.|7j it can be seen that the linear surfactants tend 
to form larger micelles on average. A summary of micelle 
statistics is given in Tab. [vj in all cases we used r^[^ — 1 [A] 
(cf Fig.|6}. 

Wetting with surfactants. Also for this train of simulations, 
the timestep size is reduced to 10 [fs]. In all other respects, 
the settings used for cylindrical droplets of pure water are 
maintained. The effect of the presence of surfactants is to 
lower the contact angle and this is possible if the tail beads, 
displaying hydrophobic behavior, interact more strongly with 
the CG graphene than the head beads P4. For this reason, it 
is assumed that the interaction strength is 2 x ecw = 0.115 
[kcal/mol]. This approach is chemically consistent since, in 
graphene, carbon atoms have no polarity. The parameter a 
for the LJ interactions with graphene beads is again 6.24 [A]. 



The results of Tab. VI indicate that the contact angle is almost 
insensitive to significant variations of surfactant concentration 
using equilibrated solutions. Since in general it is observed a 
stronger micellization process, we repeat the previous sim- 
ulations with cylindrical droplets of pure water and the sur- 
factants arranged between water and graphene as illustrated 
in Fig. [8] In this way, the solution is not yet properly equi- 
librated, but more surfactant molecules are closer to the ex- 
pected configuration, that is, near the contact line Ii22ll27ll28l . 
At 1 wt%, all the systems still do not experience any influ- 



some statistics for micelles. The critical micelle concentration 
(CMC) is the concentration at which micelles start forming. 
We calculate the CMC from the average number of free sur- 
factants lil4i . In this way, it is understood that the removal of 



ence of the presence of surfactants (see Tab. VII 1. It is worth 
noting that, in similar studies ]23[30l, low surfactant concen- 
trations resulted in a transient increase of the surface tension 
for the MARTINI model. On the other hand, the concentration 
around 4 wt% provides a basis of comparison among the dif- 
ferent systems. The results for linear surfactants prove that, in 
first approximation, the dominant factors promoting wetting 
are the length and apolarity of the hydrophobic tail. Indeed, 
comparison between surfactants L2 and L3 indicate that one 
single Ci bead is sufficient to counteract the effect of ten P4 
beads (their contact angles differ by a few degrees). We want 
to highlight the preferred arrangement of surfactant molecules 
inside the droplets. To this end, it is useful to look at the spa- 
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Figure 4: RDFs relative to the tail beads C\ of the surfactant C8E4-NaNaP4 from isobaric-isothermal dynamics. Average is taken over the 
last 100 [ns] of the whole evolution. 



Mapping 


r[K] 


Volume [A^l 


C[M] 


Average 


Variance 


Size range 


no. monomers 


CMC [M] 




298 


1'845'116 


0.158 


17.5 


81.8 


5-38 


0.0 




C8E4-NaNaP4 


305 


1'846'956 


0.158 


25.0 


244.9 


7-53 


0.0 






310 


1'847'152 


0.158 


25.0 


150.1 


8-42 


0.0 






315 


1'847'130 


0.158 


25.0 


305.0 


11-63 


0.0 




C8E4-P5 


350 


1'879'769 


0.154 


6.0 


21.1 


1-18 


5.0 


4.4 X 10-3 



Table IV: Micelle size statistics over the last 100 [ns] for the simulations of surfactant C8E4. With volume we intend the average value of the 
simulation domain; its standard deviation is always below 0.1%. The symbol C designates the concentration. 
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Figure 5: Evolution of the number of free surfactants for the C8E4-P5 
solution at 350 [K] (see Tab.llVl. 



tial distributions Py[y — ycM) and fj.(z) for Ci beads along the 
y and z axes, as the notation suggests. Here ycM is the coor- 
dinate of the center of mass of CG water. As evidenced by 



Fig. [9] linear surfactant molecules are principally located at 
the solid-liquid interface. Their organization is not uniform. 



as indicated by the statistics for micelles of Tab. VII The weak 
tendency to form layers, together with the fact that the distri- 
bution Py is strongly peaked around the contact line, give fur- 
ther support to the hypothesis that the number of hydrophobic 
units is the dominant factor for enhanced wetting. Regard- 



ing the T-shaped topology, from Fig. 10 it is seen that less 
surfactant molecules accumulate along the contact line. The 
fact that the contact angle remains almost unchanged, as com- 



pared to the linear counterparts (see Tab. VII 1, let us conclude 
that their action is more effective. In macroscopic systems, 
they should thus wet better than linear surfactants, since the 



plateau in Fig. 10 is already almost at the same height of the 
peaks near the contact line. Interestingly, a longer hydropho- 
bic tail does not lead to a lower contact angle. It thus arises 
that for T-shaped surfactants a shorter hydrophilic headgroup 
can enhance wetting more than a longer tail. The plots of 



Fig. 1 1 shows that one advantage of the T-shaped topology is 
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Figure 6: Ci-C] RDFs for the boxes with linear surfactants at 1 wt%; Left: surfactant L2; Right: surfactant L3. Similar results are obtained 
for the other equilibrated boxes with higher concentrations of surfactant molecules (cf. Tab.|Vj. 



Surfactant 


Mapping 


Volume [A^l 


no. surfactants 


C [M] - [wt%] 


Average 


Variance 


Size range 


no. monomers 


CMC [M] 


L2 


(Ci)3(P4)20 


5'825'561 
6'503'837 


34 
288 


9.7 X 10"^ - 1.0% 
7.4 X 10^2 .8 4% 


6.2 
11.5 


26.0 
54.9 


1-14 
1-85 


1.7 
1.0 


4.8 X lO"'* 
2.6 X lO"'' 


L3 


(Ci)4(P4)30 


5' 829' 183 
6'386'342 


24 
194 


6.8 X 10-^ - 1.0% 
5.0 X 10^2-8.3% 


7.4 

15.3 


6.9 
133.7 


2-11 
2-70 


0.0 
0.0 




Tl 


(Ci)3(P4)lO 


6'506'929 


520 


1.3 X 10"' -8.9% 


11.1 


54.8 


1-50 


3.9 


1.0 X lO""* 


T2 


(Ci)3(P4)20 


6'487'350 


288 


7.4 X 10-2 - 8.4% 


9.3 


41.7 


1-39 


4.2 


1.1 X lO""* 



Table V: Basic statistics regarding micelle size for the equilibration of CG water with surfactants at different concentrations. In all cases, the 
standard deviation of the volume is always around 1% of the average value. C stands for concentration. 
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Figure 7: Micelle size distributions as resulting from the simulations of equilibration for two long-chain, linear surfactants and the T-shaped 
surfactants at 8 wt% concentration. Top: For the L2 surfactant the distribution is peaked around the size of 10, while for the longer surfactant 
L3 the more frequent micelle is of size 3. Bottom: T-shaped surfactants have less micelle aggregates composed of fifteen or more surfactants. 
This occurrence has probability 0.26 for the surfactant Tl and 0. 18 for the surfactant T2. The same probabilities in the linear case are 0.33 and 
0.37 for the surfactants L2 and L3, respectively. 
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Surfactant 


Mapping 


no. water beads 


no. surfactants 


C [wt%] 


■^micelle 


Size range 


no. monomers 


d [molecules/A ] 


en 


L2 


(Ci)3(P4)20 


68'727 


43 


0.9% 


17.0 


1-14 


12.1 


0.0074 


131.1° 




















60'351 


359 


8.7% 


108.3 


1-52 


55.5 


0.0065 


130.7° 


L3 


(Ci)4(P4)30 


68'690 


28 


0.9% 


10.3 


1-10 


6.6 


0.0074 


131.1° 
















131.2° 




63'555 


239 


8.1% 


36.6 


1-66 


11.3 


0.0068 



Table VI: Characteristics of the cylindrical droplets of initial radius 130 [A] containing linear surfactants. C indicates the concentration. The 
average number of micelles is designated by A'^icelle- The cutoff distance for the extraction of micelle clusters is 7 [A] (see Fig. pi. 



Surfactant 


Mapping 


no. surfactants 


C [wt%] 


-^micelle 


Size range 


no. monomers 


d [molecules/A-'] 


en 






85 


1.1% 


41.8 


1-12 


21.6 


0.0076 


129.6° 


T 1 

Li 


(C'l)3(P4)lO 


424 


5.2% 


112.4 


1-37 


46.3 


0.0076 


118.2° 






847 


10.5% 


112.3 


1-165 


40.6 


0.0076 


114.5° 






39 


0.8% 


27.4 


1-9 


19.4 


0.0076 


130.1° 


L2 


(Ci)3(P4)20 


196 


4.2% 


75.1 


1-21 


35.4 


0.0075 


123.5° 






393 


8.4% 


91.1 


1-40 


33.9 


0.0074 


111.7° 






26 


0.8% 


15.1 


1-7 


8.5 


0.0076 


130.5° 


L3 


(Ci)4(P4)30 


132 


4.1% 


32.1 


1-21 


6.7 


0.0075 


125.0° 






265 


8.3% 


34.3 


1-43 


3.7 


0.0072 


114.9° 






530 


16.7% 


65.0 


1-206 


28.9 


0.0065 


124.7° 


Tl 


(Cl)3(P4)lO 


85 
424 


1.1% 
5.2% 


44.4 
119.1 


1-13 
1-33 


24.1 
45.4 


0.0076 
0.0076 


129.2° 
116.0° 






39 


0.8% 


28.3 


1-7 


21.0 


0.0076 


129.8° 


T2 


(Cl)3(P4)20 


196 


4.2% 


66.3 


1-20 


23.4 


0.0075 


122.1° 






393 


8.4% 


88.3 


1-36 


22.2 


0.0074 


108.3° 


T3 


(Cl)5(P4)20 


212 


5.0% 


43.8 


1-30 


9.6 


0.0075 


125.1° 



Table VII: Results for the cylindrical droplets of initial radius 130 [A] with the surfactants arranged in the initial configuration in the proximity 
of contact line; all droplets contain 68'255 water beads. C stands for concentration; A'micelle is the average number of micelles, and size range is 
the aggregation number variation. Micelle clusters are extracted using a cutoff distance of 7 [A] for all surfactants (see Fig.|6|. The simulations 
for the systems containing linear surfactants with a concentration higher than 1 wt% lasted for 100 [ns] and the analysis is performed on data 
collected under the same conditions over the last 40 [ns]. This choice is dictated by the fact that longer evolutions are necessary in order to 
reach a state near the equilibrium in these cases. The contact angle for the surfactant L3 at 16.7 wt% is not in line with the results for the lower 
concentrations suggesting that the system is not yet well equilibrated. For these longer dynamics, the stratification of water did not occur. 



faster spreading. After 25 [ns] of dynamics, the contact angles 
for T-shaped surfactants vary at most of 1.7° (surfactant T2 at 
4.2 wt%). In the same evolution period, with the exception 
of the surfactant L2 at 4.2 wt%, the variations of the contact 
angle for linear surfactants are more marked: at least of 3.9° 
(surfactant LI at 5.2 wt%) and at most of 19.2° (surfactant 
L2 at 8.4 wt%). The influence of the topology for enhanced 
spreading was already recognized in the literature l,22.,27.,28l . 



V. CONCLUSIONS 

Small concentrations of surfactants around 1 wt% can cause 
significant reductions of the contact angle. As an example, for 
the surfactant Triton® X-100, the surface tension of water of 
72.5 [mN/m] decreases linearly with concentration up to 0.03 
wt%, when the minimum of 31 [mN/m] is attained; for higher 
concentrations the surface tension remains of course constant. 
It turns out that the proper treatment of wetting phenomena 



by aqueous surfactant solutions requires conditions of difficult 
realization by the present computational capabilities. Coarse- 
grained models allowed us to move one step further toward 
a more realistic representation of such systems in terms of 
surfactant concentrations. At the typical experimental con- 
centration of 1 wt%, the coarse-grained MARTINI force field 
O'Sl marks no difference between the surfactants investigated 
here. For higher concentrations there appears that it is possi- 
ble to discriminate the wetting behavior of the various surfac- 
tants. The main instrumental conclusion of our study is that 
the length and apolarity of the hydrophobic tail determine to a 
larger extent the wetting behavior for the linear topology. In- 
stead, the length of the hydrophilic headgroup appears to be 
more relevant for the T-shaped topology. In the framework 
of our simulations, the T-shaped topology does not lead to a 
substantial decrease of contact angle vis-a-vis linear surfac- 
tants. Lennard- Jones forces are short range and the reduction 
of the contact angle is essentially driven by the accumulation 
of surfactants along the contact line. Micelle statistics demon- 
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Figure 8: (Color online) From Left to Right: Droplet with surfactant L3 at 8 wt% (see Tab. [VII^ in the initial configuration, final configuration 
without water beads and final configuration including water beads. At equilibrium, the surfactant are oriented with their head toward the fluid. 
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Figure 9: Distributions of C\ beads in space for the cylindrical droplets containing surfactants L2 at concentration 8 wt%. Top: Results 
for the droplet prepared as shown in Fig. [8] with the surfactants in the proximity of contact line. IVIany surfactants are ne ar th e contact line, 
i.e. I y — vcM l~ 90 [A], and at the solid-liquid interface. Bottom: Results obtained from an equilibrated solution (cf. Tab. VI i. The position 



and height of the peaks indicate that a significant number of surfactants is in the bulk of the droplet. Similar conclusions can be inferred from 
the outcomes for the other case studies. 



strates that linear surfactants pack more tightly, compensating 
the fact that the hydrophobic tail beads of T-shaped surfac- 
tants are on average closer to the graphitic substrate. It thus 
follows that the stronger self-assembly behavior of surfactants 
can also favor wetting. On the other hand, the weaker micel- 
lization of T-shaped surfactants can result in faster spreading. 
Finally, our parameterization of surfactants while very simple 
can encompass and provide benchmarks for several commer- 
cially available surfactants. More comparative work remains 
to be done and our results offer some insight on the complex 



interplay between micellization, spreading and wetting. 
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Figure 10: Spatial distributions of Ci tail beads for the cylindrical droplet containing surfactant T2 at 8 wt% (see Tab. |VII[ l. In this case, the 
surfactants are no longer neatly localized along the contact line. It also appears that the adsorption layer is more uniform (cf. Fig.|9|. The plots 
for the other T-shaped surfactants lead to the same conclusions. 
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Figure 1 1 : Contact angle 6 in the course of time for the cylindrical droplets of Tab. |VII| In general, T-shaped surfactants spread faster. The 
analysis uses the data collected in time intervals of 5 [ns]. 
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